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A model for J/\& suppression at a high energy heavy ion colhsion is presented. The main 
(and the only) reason for the suppression is J/\l/ inelastic scattering within hadron matter. 
The hadron matter is in the form of a multi-component non-interacting gas. Also the evolution 
of the gas, both longitudinal and transverse, is taken into account. It is shown that under 
such circumstances and with J/^ disintegration in nuclear matter added, J/'i/ suppression 
?C ' evaluated agrees well with NA38 and NA50 data. 

^ ^ PACS: 14.40. Lb, 24.10Pa, 25.75.-q 
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I. INTRODUCTION 



On the 10th of February 2000 CERN announced the discovery of the novel state of matter — the so-called 
^ . quark-gluon plasma (QGP). The existence of the QGP has been predicted upon lattice QCD calculations (for 

a review see 0| and references therein) and the critical temperature Tc for the ordinary hadron matter-QGP 
phase transition has been estimated in the range of 150 - 270 MeV (this corresponds to the broad range of the 
r^ ' critical energy density Cc ~ 0.26 — 5.5 GeV/fm^^). The upper limit belongs to a pure SU{3) theory, whereas 

f^ I adding quarks causes lowering of Tc even to about 150 MeV (cc — 0.26 GeV/fm~^ respectively). Since NA50 

Collaboration estimates for the energy density obtained in the central rapidity region (CRR) give the value of 
3.5 GeV/fm~^ for the most central Pb-Pb data point g], it is argued that the region of the existence of the 
Pj . QGP has been reached in Pb-Pb collisions at CERN SPS. 

^ i' Soon after the announcement, the paper entitled "Evidence for deconfinement of quarks and gluons from the 

J/^ suppression pattern measured in Pb-Pb collisions at the CERN-SPS" appeared 0]. As the title suggests, 
the main argument for the QGP creation during Pb-Pb collisions at the CERN SPS is the observation of the 
r^ ', suppression of J/^ relative yield. J/* suppression as a signal for the QGP formation was originally proposed 

by Matsui and Satz W. The clue point of the NA50 Collaboration paper is the figure (denoted as Fig. 6 



there), where experimental data for Pb-Pb collision values of ^^ ''^''' (the ratio of the J/^ to the Drell-Yan 
production cross-section times the branching ratio of the J/^ into a muon pair) are presented together with 
^ . some conventional predictions. Here, "conventional" means that J/^ suppression is due to J/vP absorption in 

ordinary hadron matter. Since all those conventional curves saturate at high transverse energy Et, but the 
experimental data fall from Et — 90 GeV much lower and this behaviour could be reproduced on the base of 
J/^ disintegration in QGP H, it is argued that QGP had to appear during most central Pb-Pb collisions. In 
general, the immediate reservation about such reasoning is that besides those already known (see e.g. [|6|-[l0[ and 
references [12-15] in B), there could be a huge number of different conventional models, so until this subject is 
cleared up completely, no one is legitimate to ruled out the absorption picture of J/^ suppression. 

In the following talk, the more systematic and general description of J/^ absorption in the framework of 
statistical analysis will be presented. The main features of the model are [Q: 

• 1. a multi-component non-interacting hadron gas appears in the CRR instead of the QGP; 

• 2. the gas expands longitudinally and transversely; 

• 3. J/^ suppression is the result of inelastic scattering on constituents of the gas and on nuclcons of 
colliding ions. 



*Talk based on the work done in collaboration with L.TurkotJ. 

^Presented at the XLI Cracow School of Theoretical Physics "Fundamental Interactions", Zakopane, Poland, June 2-11 
2001. 
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II. THE TIMETABLE OF EVENTS IN THE CRR 

For a given A-B collision i = is fixed at the moment of the maximal overlap of the nuclei (for more details 
see e.g. P]). As nuclei pass each other charmonium states are produced as the result of gluon fusion. After half 
of the time the nuclei need to cross each other (t ~ 0.5 fm), matter appears in the CRR. It is assumed that 
the matter thermalizes almost immediately and the moment of thermalization, to, is estimated at about 1 fm 
P,0|. Then the matter begins its expansion and cooling and after reaching the freeze-out temperature, Tfo., it 
ceases as a thermodynamical system. The moment when the temperature has decreased to Tfo. is denoted as 
tf,o.- Since the matter under consideration is the hadron gas, any phase transition does not take place during 
cooling here. 

For the description of the evolution of the matter, relativistic hydrodynamic is explored. The longitudinal 
component of the solution of hydrodynamic equations (the exact analytic solution for an (l+l)-dimensional 
case) reads (for details see e.g. |0,n2]) 



t •, ^oTo , . n%TQ z 

T T t 



where r = \Jt'^ — z^ is a local proper time, v^ is the z-component of fluid velocity [z is a collision axis) and sq 
and rP^ are initial densities of the entropy and the baryon number respectively. For ub = and the uniform 
initial temperature distribution with the sharp edge at the border established by nuclei radii, the full solution 
of (3+l)-dimensional hydrodynamic equations is known p3[ . The evolution derived is the decomposition of the 
longitudinal expansion inside a slice bordered by the front of the rarefaction wave and the transverse expansion 
which is superimposed outside of the wave. Because small but nevertheless non-zero baryon number densities 
are considered here, the above-mentioned description of the evolution has to be treated as an assumption in 
the presented model. The rarefaction wave moves radially inward with a sound velocity Cj. The sound velocity 
squared is given by c^ = ^ and can be evaluated numerically ||lj,|l5[ . 

HI. THE MULTI-COMPONENT HADRON GAS 

For an ideal hadron gas in thermal and chemical equilibrium, which consists of I species of particles, energy 
density e, baryon number density n^, strangeness density ng and entropy density s read (h = c= \ always) 
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J2{2s, + 1) dp ^^-^ , (2c) 

' |-oo 4 (£;^_^^)exp|^^^| 

^(2..-fl)/ dpP- I -I, (2d) 

■'' ^' (exp{^V^}+g,) 



where Ei = {mj -\- p^Y'"^ and m^, Bi, Si, pi, Si and gi are the mass, baryon number, strangeness, chemical 
potential, spin and a statistical factor of specie i respectively (an antiparticle is treated as a different specie) . 
And Pi = BipB -f Sips, where fiB and fis s^re overall baryon number and strangeness chemical potentials 
respectively. 

To obtain the time dependence of temperature and baryon number and strangeness chemical potentials one 



has to solve numerically equations (2b - 2d) with s, ub and ns given as time dependent quantities. For s{t) and 
nB{T) expressions (|l|) are taken and ns — since the overall strangeness equals zero during all the evolution 
(for more details see M]). 
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IV. J/* ABSORPTION IN THE EXPANDING HADRONIC GAS 

As it has been already mentioned in Sect. |ll| charmonium states are produced in the beginning of the coUision, 
when nuclei overlap. For simplicity, it is assumed that the production of cc states takes place at t = 0. To 
describe J/\l/ absorption quantitatively, the idea of Q is generalized to the case of the multi-component massive 
gas, here. Since in the CRR longitudinal momenta of particles are much lower than transverse ones (in the 
c.m.s. frame of nuclei), J/^ longitudinal momentum is put at zero. Additionally, only the plane z = is under 
consideration, here. For the simplicity of the model, it is assumed that all charmonium states are completely 
formed and can be absorbed by the constituents of a surrounding medium from the moment of creation. The 
absorption is the result of a cc state inelastic scattering on the constituents of the hadron gas through interactions 
of the type 



D + D + X, 



(3) 



where h denotes a hadron, I? is a charm meson and X means a particle which is necessary to conserve the 
charge, baryon number or strangeness. 

Since the most crucial for the problem of QGP existence are Pb-Pb collisions (see remarks in Sect. B), the 
further considerations are done for this case. So, for the Pb-Pb collision at impact parameter b, the situation 
in the plane z = is presented in Fig.^, where Seff means the area of the overlap of the colliding nuclei. 




FIG. 1. View of a Pb-Pb collision at impact parameter b in the transverse plane (z 
overlap is hatched and denoted by Seff- 



0). The region where the nuclei 



Additionally, it is assumed that the hadron gas, which appears in the space between the nuclei after they 
have crossed each other, also has the shape of Seff at to in the plane z = 0. Then, the transverse expansion 
starts in the form of the rarefaction wave moving inward S^ff at ip- 

From the considerations based on the relativistic kinetic equation (for details see |ll|,||), the survival fraction 
of J/\I/ in the hadron gas as a function of the initial energy density eg in the CRR is obtained: 



/ 



J^h.gXeo) = / dpT g{pT, eo) • exp ■ 



dt 
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where the sum is over all taken species of scatters (hadrons), p'^ = {E,pt) and gf ~ {E[, q) are four momenta 
of J/^ and hadron specie i respectively, Ui states for the absorption cross-section of J /^ — hi scattering, Vrei,i 
is the relative velocity of hi hadron with respect to J/^ and M and rrn denote J/^ and hi masses respectively 
(M = 3097 MeV). The function g{pT,£o) is the J/^ initial momentum distribution. It has a gaussian form 
and reflects gluon multiple elastic scatterings on nucleons before their fusion into a J/^> in the first stage of 
the collision [|l6|-[19|. The upper limit of the integration over time in (0), namely t final is the minimal value 
of (tesc) and tfo,. The quantity (tesc) is the average time of the escape of J/^'s from the hadron medium for 
given values of b and J/'i velocity v = pr/E. Note that the average is taken with the weight 



Pj/^ir) 
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(5) 
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where Tj^(s) — J dzpj^{s, z) and pa{s,z) is the nuclear matter density distribution. For the last quantity, the 
Woods-Saxon nuclear matter density distribution with parameters from ||20(| is taken. In the integration over 
hadron momentum in (0) the threshold for the reaction (Bl) is included, i.e. ai equals zero for (p"^ + gj')^ < 
(2tou + mx)"^ and is constant elsewhere (mu is a charm meson mass, tod — 1867 MeV). Also the usual 
Bose-Einstein or Fermi-Dirac distribution for hadron specie i is used in M) 

fMt)-Mci,t)^ ':;+^| — . (6) 

exp|^T(t}— /+ft 

As far as ai is concerned, there are no data for every particular J/^ — hi scattering. Therefore, only two types 
of the cross-section, the first, (Tj,, for J/^ — baryon scattering and the second, am, for J/^ — meson scattering 
are assumed, here. In the model presented at is put at (ij/^iq — J/^ — Nucleon absorption cross-section. 
And values of (Jj/tpn in the range of 3-5 mb are taken from p-A data pl|-p3[. The J/^ — meson cross-section 
CTm is assumed to be 2/3 of ct;,, which is due to the quark counting. 

As it has been already suggested ^^ also J/^ scattering in nuclear matter should be included in any J/^ 
absorption model. This could be done with the introduction of J/^ survival factor in nuclear matter |23-|29] 

Mi.m.(eo) = exp{-crj/^jv/9oi} , (7) 

where po is the nuclear matter density and L the mean path length of the J/\I' through the colliding nuclei. 
The length L is expressed by the following formula pfl]: 



P^m = -l-Jd'sTA{s)TBis- b) 



±^TA{s) + ^Ts{s-b) 



(8) 



where TAsib) = J d^sTAisjTsis- b). 

Since J/vP absorptions: in nuclear matter and in the hadron gas, are separate in time, J/^ survival factor 
for a heavy-ion collision with the initial energy density eo, could be defined as 

A/'(eo) = Mn.rn. (co) ' J^h.g. (eo) ■ (9) 

Note that since right sides of (0) and (|7|) include parts which depend on impact parameter b and the left sides 
are functions of eq only, the expression converting the first quantity to the second (or reverse) should be given. 
This is done with the use of the dependence of eq on the transverse energy Et extracted from NA50 data M 
(for details see Q). 

To make the model as much realistic as possible, one should keep in mind that only about 60% of J/^'s 
measured are directly produced during collision. The rest is the result of x (^ 30%) and -0' (~ 10%) decay [ p7[ . 
Therefore the realistic J/^ survival factor could be expressed as 

AA(eo) = 0.6Uj/4eo) + 0.3AAx(eo) + 0.1A/'^.(eo) , (10) 

where Afj/^{eo), J\f^{eo) and A/"^' (eg) are given also by formulae (^, (0) and (^) but with g{pT,eo), ctj/^n and 
M changed appropriately (for details see y). 

To complete, also values of cross-sections for x~baryon and ip' — baryon scatterings are needed. For simplicity, 
it is assumed that both these cross-sections are equal to J/^ — baryon one. Since J/^ is smaller than % or ?/;', 
this means that J/^ suppression evaluated according to (HOh is underestimated here. 

V. RESULTS 

To calculate formula (0) initial values sq and n^ are needed. For the estimation of the last quantity experi- 
mental results for S-S ||28] and Au-Au |2^j32l collisions are explored. The S-S data give n^ = 0.25 fm~'^ and 
Au-Au data n^ = 0.65 fm"'^. These values are for central collisions, for more peripheral ones the initial baryon 
number density should be lower. To simplify numerical calculations, n^ is kept constant over the whole range 
of b. But, to check the possible dependence on n^, also the much lower value n^ = 0.05 fm~'^ is examined. To 
find So J (|2^ - |2c|) with respect to T, /ig and fiB are solved. Then, having put Tq, /ig and fi% into (pq), sq is 
obtained. The last parameter of the model is the freeze-out temperature. Two values Tfo. — 100, 140 MeV are 
taken here and they agree well with estimates based on hadron yields [p3 . 

Firstly, formula ( [lO| ) with the use of (^) only instead of (g) is calculated for the case of (H-l)-dimensional 
expansion and results are presented in Figs.0-pl 
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FIG. 2. J/^ suppression in the longitudinally 
expanding hadron gas for n% — 0.25 fm~^ and 
Tf.o. = 140 MeV: (a) n = 3 mb, a^ = 2 mb; (b) 
(7i, = 4 mb, am = 2.66 mb; (c) at = 5 mb, cr™, = 3.33 
mb; (d) CT(, = 6 mb, (Tm = 4 mb . The black squares 
correspond to the NA38 S-U data [Bl|, the black tri- 
angles correspond to the 1996 NA50 Pb-Pb data, the 
white squares to the 1996 analysis with minimum bias 
and the black points to the 1998 analysis with mini- 
mum bias fsl. 




£o [GeV/fm^] 
FIG. 3. Same as Fig.H (except case (d), which is 
not presented here) but for T/.o. = 100 MeV. 



For comparison, also the experimental data are shown in Figs. 



The experimental survival factor is defined 



as 



•J^ exv 



^ exp 



DY 



J/i' 



(11) 



where 



AB(pp) 
DY 



is the ratio of the J/^ to the Drell-Yan production cross-section in A-B(p-p) interactions 



times the branching ratio of the J/^ into a muon pair. The values of the ratio for p-p, S-U and Pb-Pb are taken 
from 1^,^,^. The results of numerical evaluations of (|lO|) agree with the data well qualitatively for greater 
baryonic cross-section a^, and (or) for the lower freeze-out temperature Tf^., except the last point measured for 
the highest Et bin. 

Now the full (3-M)-dimcnsional hydrodynamic evolution (in the form described in Sect. O) and J/^ absorption 
in nuclear matter will be taken into account and results are depicted in Figs.Q-||. To draw also S-U data together 
with Pb-Pb ones, instead of multiplying Mh.g. by Mn.r 
hadron gas survival factor" is defined as 



Mexp is divided by A/",!.™., i.e. "the experimental J/^ 



Mex-p = CXp {cr,//^Ar/Ooi} ' A/'e: 



(12) 



and values of this factor are drawn in Figs.0-ra as the experimental data. The quantity tq is a parameter from 
the expression for a nucleus radius Ra = tq • As. The value of Ka is used to fix the initial position of the 
rarefaction wave in the evaluation of (tesc) in ^^ ^f Koo f,,,-Tio^ ^i,f +kQ+ +Ko n^aa ^f ^o — n n^; f,Tn-3 
differ substantially from that of n^ — 0.25 fm 



|). It has turned out that the case of n^ — 0.05 fm ■^ does not 
so curves for n^ = 0.05 fm~^ are not depicted in Figs.Era. 
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0.0 0.5 1.0 1.5 2.0 2.5 
Eo [GeV/frn] 
FIG. 5. Same as Fig.^ but for a^ -- 
cTm = 3.33 mb. 



5 mb and 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
£o [GeV/fm^] 
FIG. 4. J/^ suppression in the longitudinally 
and transversely expanding hadron gas for the 
Woods-Saxon nuclear matter density distribution and 
0-^=4 mb, am = 2.66 mb and Tf,o. = 140 MeV. The 
curves correspond to n% = 0.25 fm""^, Cg — 0.45, 
ro = 1.2 fm (solid), n^ = 0.65 fm-^ Cs = 0.46, 
ro = 1.2 fm (dashed) and n% = 0.25 fm-^ Cs = 0.45, 
ro = 1.12 fm (short-dashedV The black squares rep- 
resent the NA38 S-U data pn , the black triangles rep- 
resent the 1996 NA50 Pb-Pb data, the white squares 
the 1996 analysis with minimum bias and the black 
points the 1998 analysis with minimum bias [H], but 
the data are "cleaned out" from the contribution of 
J/^ scattering in the nuclear matter in accordance 
with (^. 

Having compared Figs. 0-0 with Figs.p[«, one can see that adding the transverse expansion changes the final 
(theoretical) pattern of J/^ suppression qualitatively. Now the curves for the case including the transverse 
expansion are not convex, in opposite to the case with the longitudinal expansion only, where the curves are. 
As far as the pattern of suppression is concerned, theoretical curves do not fall steep enough at high eg to 
cover the data area. But for some choice of parameters, namely for at somewhere between 4 and 5 mb and 
for ro — 1.12 fm, a quite satisfactory curve could have been obtained. Precisely, again only the highest Et 
bin point falls down of the range of theoretical estimates completely. But the error bar of this point is very 
wide. And also the contradiction in positions of the last three points of the 1998 data can be seen. This means 
that the high Et region should be measured once more with the better accuracy to state definitely whether the 
abrupt fall of the experimental survival factor takes place or not there. 

To support the conclusion, main results from Figs.^-^ are presented in Fig.|^ again. The original data |^] for 

g ^PbPb g g.PP 

'^^pb'pb and J/^ survival factors given by (flfl) multiplied by ''^p'^'*' and now as functions of Et are depicted 
there (for details see ffl). 
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FIG. 6. J/*!/ survival factor times 



-—m> 



in the longitudinally and transversely expanding hadron gas for the 

' DY 

Woods-Saxon nuclear matter density distribution and n^ = 0.25 fm"'', Tf,o. = 140 MeV, Cs — 0.45 and ro = 1.2 fm. The 
curves correspond to (Jj, = 4 mb (solid) and a;, = 5 mb (dashed). The black triangles represent the 1996 NA50 Pb-Pb 
data, the white squares the 1996 analysis with minimum bias and the black points the 1998 analysis with minimum bias 



As it has been already mentioned, the main disagreement with the data appears in the last experimental point 
of the 1998 analysis. 

The charmonium — baryon inelastic cross-section is the most crucial parameter in the model presented. So 
to be sure what is the exact result of J/^ absorption, one should know how this cross-section behaves in the 
hot hadron environment. The newest estimations of tt + J/^, p + J/^ and J/'f + N cross-sections at high 
invariant collision energies |33|j34| give values of at, and cr,„ of the same order as assumed here. Also it should 
be stressed that the charmonium — hadron inelastic cross-sections are considered as constant quantities in 
the model presented. But, as the results of just mentioned papers suggest, they are not constant at all. The 
cross-sections are growing functions of the invariant collision energy yG. Therefore, one could think naively 
that the increase of eq (or in other words Et) causes the increase of the invariant collision energy ^/s on the 
average and further the increase of the charmonium — hadron inelastic cross-sections. In this way, the line 
describing J/^ survival factor could be placed close to the sohd curve of Fig. ^ for low eo (Et), and then, as the 
charmonium — hadron inelastic cross-sections increase, this line would go closer to the dashed curve of Fig.|^ 
for high eo {Et)- So, the experimental pattern of J/^ suppression could be recovered. 

As the last remark, the author would like to call reader's attention to Fig. 5 of |^. The solid and dotted 
curves in that figure are exactly the same as curves just presented in Fig. 0. But the distortion in the QGP is 
the main reason for J/^ suppression in |35[ ! It is also stated there, that results shown "provide evidence for 
the production of the quark-gluon plasma in central high-energy Pb-Pb collisions" . 

Therefore, the final conclusion is the following: the NA50 PB-PB data provides evidence for the production 
of the thermalized and confined hadron gas in the central rapidity region of a PB-PB collision. This implies that 
J/vf* suppression is not a good signal for the quark-gluon plasma appearance at a central heavy-ion collision. 
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